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SUMMARY 

This primer provides documentation for using MSC NASTRAN In analyzing 
rotating flexible blades. The analysis of these blades Includes geometrically 
nonlinear (large displacement) analysis under centrifugal loading, and fre- 
quency and mode shape (normal modes) determination. The geometrically non- 
linear analysis using NASTRAN Solution sequence 64 is discussed along with the 
determination of frequencies and mode shapes using Solution Sequence 63. A 
sample problem with the complete NASTRAN Input data Is Included. Items unique 
to rotating blade analyses, such as setting angle and centrifugal softening 
effects are emphasized. 


1.0 INTRODUCTION 

The purpose of this primer Is to document the use of MSC NASTRAN In ana- 
lyzing rotating flexible blades. The analysis of rotating flexible blades such 
as compressor and turboprop blades, often requires complex procedures Including 
geometrically nonlinear (large displacement) analysis and frequency and mode 
shape determination. The objective in performing such analyses Includes the 
prediction of steady state deflections and stresses under centrifugal forces, 
the generation of data for constructing Campbell diagrams (plots of frequency 
versus rotational speed), and the provision of modal data for use in flutter 
calculations. In performing these analyses, and in modeling the complex geom- 
etries and material properties of the blades (fig. 1), finite element (F.E.) 
computer programs typically are used. NASTRAN is particularly well-suited 
because of Its ability to compute steady-state displacements from Its geometri- 
cally nonlinear analysis capabilities, and then, to use those results for a 
subsequent normal modes analyses. 

This primer also can be used to bridge the gap between the theoretical 
modeling tactics presented In (ref. 1) for creating the model, and the prac- 
tical application of NASTRAN for large displacement or normal modes analysis. 
Reference 1 provides a complete discussion of blade modeling strategies and 
documents the capabilities and limitations of various NASTRAN elements. The 
performance of a broad range of mesh configurations and element types also are 
evaluated for different blade related parameters such as camber, twist, sweep, 
and rotational speed. 

The computation of steady-state displacements, frequencies, and modes 
shapes of flexible rotating blades requires that two NASTRAN Solution sequences 



be run. First, a large displacement analysis Is run using NASTRAN Solution 
Sequence 64. This solution sequence performs large displacement analysis on 
the rotating blade, computes steady-state displacements and stresses (due to 
the rotational effects), and then stores the blade's final stiffness and mass 
matrix In a database. Experience has shown (ref. 2) that a large displacement 
analysis Is required because the blades are relatively flexible and normally 
deflect considerably under the centrifugal forces (fig. 2). 

Following the large displacement analysis, the frequencies and mode shapes 
are computed using Solution Sequence 63. (A typical plot demonstrating the 
variation In natural frequencies with rotational speed Is shown In fig. 3.) 

This solution sequence computes the modal parameters from the final mass and 
stiffness matrices which were computed during the Solution 64 run. The struc- 
tural property matrices that correspond to the blade in its deformed position 
must be used so that the effects of centrifugal stiffening, and other elements 
which will be discussed, are Included In the normal modes analysis. 


2.0 DESCRIPTION OF SAMPLE PROBLEM 

The sample problem In figure 4 Is provided In order to demonstrate the 
procedures required for large displacement and normal modes analyses. All of 
the requirements for these analyses are exhibited In this sample problem even 
though this problem does not have the mesh complexity required of typical flex- 
ible blades. The mesh used In this sample problem Is for demonstration pur- 
poses and Is over simplified compared to typical analyses. Although there are 
no aerodynamic loads included in this problem, they can be Included, and would 
be combined with the centrifugal loads for computing the steady-state 
displacements . 

The sample problem consists of a rotating, swept back, flat plate. The 
plate Is modeled with three plate elements (CQUAD4) connected by 8 grid points. 
The grid points are defined In a local rectangular coordinate system, which is 
rotated 30° from the axis of rotation In the X/Y plane (see section 6.0 for 
a discussion on the C0RD2R card). It should be noted that the selection of 
the orientation of axes needs to be consistent with the directions that are 
specified for the centrifugal softening terms and the rotational speed 
(section 3.0). 


3.0 GEOMETRICALLY NONLINEAR ANALYSIS/MSC NASTRAN SOLUTION 64 

This section provides a discussion regarding the construction of the 
Solution 64 Input data deck. The sample data deck Is comprised of four compo- 
nents; the CRAY Job Control Language, the NASTRAN Executive, the Case Control, 
and Bulk Data, and Is Included In appendix A. The following sections discuss 
step by step the data cards Included In the appendix. 


3.1 CRAY Job Control Language 

The CRAY Job Control Language for submitting jobs to the CRAY (at NASA 
Lewis Research Center) and running NASTRAN Is given at the beginning of the 
data deck. The amount of time and memory Indicated on the "JOB" card Is based 
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on the number of degrees of freedom used for the model and the number of Itera- 
tions (section 3.3) specified for the nonlinear analysis. Faster turnaround 
time on the computer Is accomplished when these allocations are minimized. 

On the "NASTRAN" card (see ref. 3 for a complete description), a temporary 
database named “BLADE" Is specified. This database Is used for storing the 
mass and stiffness matrices for subsequent use in the Solution 63 normal modes 
analysis . 


3.2 Solution 64 Executive 


Solution 64 Is well suited to the large steady-state displacement analysis 
of rotating structures except for the fact that corlolls forces and centrifugal 
softening terms (see appendix C for a derivation of centrifugal softening) are 
not automatically Included. Since the corlolls forces are velocity related, 
they do not have any Influence on the large displacement analysis and do not 
need to be accounted for. Furthermore, from previous studies It has been 
determined that corlolls forces have negligible effects on thin, flexible blade 
frequencies (ref. 4), so they normally do not need to be Included in the Solu- 
tion 63 analysis. However, the centrifugal softening terms do need to be 
Included In the analysis and must be added via DMAP (Direct Matrix) 
programming. 


The centrifugal softening terms (-w 2 M) are added to the global stiffness 
matrix for all grid points for the translational degrees of freedom In the 
direction of the two axes perpendicular to the axis of rotation. The softening 
terms are input Into the Solution 64 analysis using the NASTRAN DMAP included 
with the sample problem given In appendix A. In DMAP Alter 77, the -w 2 
coefficient first Is computed from the rotational speed specified on the RPM 
parameter card which is provided In the bulk data deck (see section 3.4). 

This coefficient then is added Into the matrix "K1GM". All terms In this 
matrix are zero except for the diagonal elements corresponding to the two 
translational degrees of freedom perpendicular to the axis of rotation. These 
elements are set equal to -w 2 . The matrix "Kl" which specifies the two per- 
pendicular axes is defined using DMI cards In the bulk data deck (section 3.4). 
The K10M matrix Is multiplied by the global mass matrix to form the -w 2 M 
softening matrix and then the resulting matrix, KSQFT, is stored In a database. 
In alter 188 KSOFT Is retrieved from the database ( "DBFETCH. . . ") , and Is 
added to the global stiffness matrix (K3J) In every Iteration that takes place 
In the large displacement analysis. 


3.3 Solution 64 Case Control 

The cards In the CASE CONTROL Deck are used to specify the problem titles, 
the type of printed output, and the number of iterations that are to be carried 
out In the large displacement analysis. 

For the sample problem the applied forces at all of the grid points 
( OLOAD=ALL) , and the resulting displacements ( D ISP-ALL) and reactions 
(SPCF0=ALL), are printed in the last three Iterations. It Is recommended that 
displacements and reactions be printed In at least the last few of Iterations 
so that convergence can be monitored. It may be desirable to print displace- 
ments In the first Iteration so that the nonlinearity of the blade's response 
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can be assessed. Acceptable convergence Is achieved when both the displacement 
changes between Iterations and the force unbalance between the applied centri- 
fugal forces and the Internal element forces are small. The force unbalance 
at each of the unconstrained grid points Is printed along with the reactions 
at the constrained points by using the command "SPCF0=ALL" . 

There are eight Iterations, or "subcases," specified for this sample 
problem. In the first subcase a linear analysis Is performed. In the second 
subcase, the displacements from the first subcase are used to form a differen- 
tial stiffness matrix which Is then used to compute a new set of displacements. 
In subsequent subcases the differential stiffness matrix is updated using the 
resultant displacements and a new set of displacements Is computed. (Details 
of the theory underlying the Iterative procedure and results for a flexible 
turboprop blade are provided In (refs. 2 and 5).) 

When actual blades having large numbers of degrees of freedom are analyzed 
the cost of running a large number of Iterations can be significant. To mini- 
mize the CPU time and cost, the number of Iterations should be kept to a mini- 
mum. The best way to optimize the number of Iterations Is to specify a minimum 
number of Iterations, and then If the solution has not converged, use NASTRAN 
"restart" (ref. 6) to resume the analysis. 


3.4 Solution 64 Bulk Data Deck 

The primary function of the bulk data deck Is to describe the blade geom- 
etry, boundary conditions, material properties, and loads. Details of blade 
modeling techniques and the bulk data cards required for the model description 
are presented In (ref. 1), and In the NASTRAN User's Manual (ref. 6). In addi- 
tion to describing the model, the bulk data deck Is used to specify the rota- 
tional speed (RFORCE and PARAM RPM), and the matrix "XI" discussed In 
section 3.2. 

NASTRAN automatically computes a centrifugal force field whenever an 
"RFORCE" card is used to specify a rotational speed. The centrifugal force 
field Is computed by using the blade's geometrical and mass properties defined 
In the bulk data deck, and the rotational speed specified on the "RFORCE" card. 
The rotational speed Is also included on a parameter card, "PARAM RPM," This 
card Is used by the DMAP Alters for computing the value -w^ for the centri- 
fugal softening terms. 

The matrix K1 Is set up using DMI Bulk Data cards. In the first DMI 
card, the matrix Is defined as being symmetric and of size 6 by 6. The size of 
this matrix Is always 6 by 6 corresponding to the number of degrees of freedom 
at a grid point. In the next two DMI cards the elements of the matrix are 
defined as zero except for the diagonal terms corresponding to the two trans- 
lational directions perpendicular to the axis of rotation which are set equal 
to 1.0. For this problem the Y and Z axis are the axes that are perpen- 
dicular to the axis of rotation. In alter 77 (section 3.2) the K1 matrix 
will be expanded so that the softening terms are applied at every grid point 
In the model. 

Several Issues which concern the blade model and are. addressed In (refs. 1 
and 7), are discussed below as they are relevant to the large displacement ana- 
lysis of flexible blades. The first Issue concerns the method of formulation 



for the element mass matrices. Most of the elements available In NASTRAN per- 
mit the user to utilize a lumped or consistent mass matrix, but since the form- 
ulation for the centrifugal softening terms, and the centrifugal force field, 

Is based on a lumped mass representation, a lumped mass matrix should also be 
used for the elements. Furthermore, no clear advantages have been found for 
utilizing a consistent mass formulation. A lumped mass matrix Is computed by 
default In NASTRAN. 

The next Issue regards the lack of stiffness In the In-plane, normal rota- 
tion for plate and shell elements. This condition can present problems when 
adjacent elements lie In the same plane (coplanar). Since the elements have 
zero computed stiffness in the normal rotation, and the elements are coplanar, 
the accumulated rotational stiffness may end up being very close to zero. When 
this occurs the stiffness matrix becomes singular and the analysis falls. To 
circumvent this problem, these "small 1 ' rotational stiffnesses can be con- 
strained with SPC's. This solution Is sensible since the normal rotational 
stiffness usually Is relatively stiff, thus It Is reasonable to fully con- 
strain the rotation. For nonlinear problems this approach Is not feasible 
because the large displacements may deflect the elements such that elements 
that start out noncoplanar become coplanar. When this happens, singularities 
that did not exist In earlier iterations arise, and the solution falls In the 
resultant Iteration. A feature In NASTRAN for overcoming this problem Is the 
K6R0T parameter. This parameter adds artificial rotational stiffness at the 
element level so that even If elements are coplanar, the global stiffness 
matrix will not have any singularities. The effect of adding this rotational 
stiffness produces results similar to when the rotation is fully constrained 
with SPC's (see ref. 1 for details of comparison). The advantage of using the 
K6R0T parameter over using SPC's Is that it Is difficult to determine before- 
hand where SPC constraints are required. 

The final Issue concerns the constraint at the base of the blade model. 

For most blade analyses performed thus far, the base of the blade has been 
fully constrained. In reference 7, It was shown that base flexibility can have 
a significant effect on steady-state displacements, frequencies, and mode 
shapes. Therefore, whenever there is information on the blade's base support 
flexibility. It should be incorporated into the blade model. 


4.0 NORMAL MODES ANALYSIS/MSC NASTRAN SOLUTION 63 

The NASTRAN Solution 63 data deck for the normal modes analysis is given 
in appendix B. This data deck uses the Identical model description that was 
used In the Solution 64 deck. The Solution 63 deck was created by duplicating 
all of data cards used in the Solution 64 deck for defining the blade model, 
removing the cards associated with the large displacement analysis, and then 
adding the necessary cards for the Solution 63 normal modes analysis. 


4.1 CRAY Job Control Language 

The amount of time and memory specified for the Solution 63 run is based 
on the number of degrees of freedom used In the analysis and the number of 
modes requested. As with the Solution 64 analysis, faster turnaround time on 
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the computer will be accomplished when the time and memory allocations are min- 
imized. It also should be noted that the "NASTRAN" command accesses the same 
database ("BLADE") In the Solution 63 run as in the Solution 64 run. 


4.2 Solution 63 Executive 

The primary function of the Executive Control Deck is to request that 
Solution Sequence 63 be run (SOL 63). 


4.3 Solution 63 Case Control 

The Case Control deck specifies that the mass and stiffness matrices are 
to be retrieved off of the database (SEKR and SEMR) , and that the method to be 
used for the elgensolutlon Is on "EIGR" card 10 In the bulk data deck 
(METH0D=10). The "DISP=ALL" card specifies that the values of the computed 
mode shapes at all of the grid points are to be printed out. 


4.4 Solution 63 Bulk Data Deck 

The model description provided In the Solution 63 bulk data deck is not 
used for creating a new mass and stiffness matrix since these matrices are 
obtained from the database. These cards are required In the Data Deck for 
bookkeeping purposes. 

The "EIGR" Card is used to specify the method for extracting the blade's 
frequencies and mode shapes. For this sample problem the Inverse power method 
is chosen. This method of solution Is recommended for flexible blades with a 
large number of degrees of freedom, and when there Is Interest only in the 
first few modes. For this blade, the frequency range of Interest Is specified 
to be between 50 and 3000 Hz with an estimate that there will be 7 modes In 
this range. The analyst must use care with these specifications because. If 
the range Is made too large, there may be more modes In the range than wa$ 
estimated and some of the modes may be missed, or the modes that are computed 
may be computed Inaccurately. Conversely, If the range Is made too small some 
of the modes may be outside of the range and not computed at all. In 
section 3.1.2 of the NASTRAN USER'S MANUAL (ref. 6) a description is given of 
the elgensolutlon summary which Is printed along with the NASTRAN output when 
the Inverse power method Is used. This summary provides the analyst with 
Information concerning the adequacy of the elgensolutlon and the reason for 
solution termination. The third reason for termination listed in section 3.1.2 
(Indicated by "All eigenvalues found In the frequency range specified") is most 
desired as It indicates that the solution was completed correctly. If this 
Is not the reason for termination, the frequency range or the estimated number 
of modes In the range specified on the EIGR card likely will need to be 
adjusted . 


5.0 NONROTATING BLADE NORMAL MODES ANALYSIS 

The data deck for the Solution 63 analysis can be used, when It is neces- 
sary to obtain the nonrotating blade frequencies and mode shapes (NASTRAN 
Solution Sequence 3 can also be used). When utilizing nonrotating blades In 
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the normal modes analysis SELA = 1 and the SEMA =1 in the Solution 63 case 
control deck need to be changed to SELA = all and SEMA - all. When these 
changes are made, NASTRAN computes a new mass and stiffness matrix Instead of 
pulling the mass and stiffness matrices off of the database which Is done when 
SELA and SEMA are set to 1. Thus, the rotational effects computed in the 
large displacement analysis are not Included In the property matrices since 
they are computed from data supplied for the nonrotating blade. 


6.0 CHANGES IN ROTATIONAL SPEED 

To obtain the steady state displacements, frequencies, and modes shapes 
when the blade is spinning at a new speed, both the large displacement and 
normal modes analyses must be rerun. The large displacement analysis (Solu- 
tion 64) needs to be rerun because the steady-state position of the blade 
changes with rotational speed. Furthermore, updated mass and stiffness 
matrices, which reflect the effects of the current rotational speed, need to 
be transferred to the normal modes analysis performed in Solution 63. The 
effect of rotational speed on blade frequencies for a typical blade was shown 
In figure 3. 

In the Solution 64 run the rotational speed is altered by changing the 
speed on the RFORCE and PARAM RPM cards. In addition to changing the blade’s 
speed, the blade's angle of attack normally has to be adjusted when the blade 
Is operating at a new rotational speed. Once the correct angle of attack Is 
determined, either the axis that the blade is rotating about must be changed 
or the blade Itself must be rotated. Due to the way the centrifugal softening 
terms are applied It Is simpler to rotate the blade than to change the axis of 
rotation. To Implement a change In the angle of attack, the entl re blade can 
be rotated by defining the blade's geometry in a new coordinate system. This 
method of rotating the blade is convenient because the coordinates of the blade 
on the "GRID" cards do not have to be changed. Instead, the coordinate system 
In which the blade is described Is rotated. Referring to figure 5, the blade 
cross section shown Is rotated an angle a by prescribing that the blades 
coordinates are In reference to the X', Y‘, Z' system Instead of from the 
original X, Y, Z system. The X 1 , Y', Z 1 system Is defined by using the 
"CORD2R" card which is described In section 7.3.2 of the NASTRAN Primer 
(ref. 8). The blade's grid points are designated to be In the primed coordi- 
nate system by changing the Identification number of the coordinate system 
using a GRDSET bulk data card. It should be noted that all other items, such 
as the blade's rotational speed, continue to be in reference to the original, 
basic coordinate system. 


7.0 COMBINED SOLUTION 64/SOLUTION 63 ANALYSIS 

A combined large displacement, frequency, and mode shape analysis can be 
performed in the Solution 64 Sequence. This analysis is completed by adding 
NASTRAN DMAP Alters to the Solution 64 data deck (see appendix D) . The func- 
tion of these DMAP alters Is to access the eigenvalue extraction routines that 
solve for the frequencies and mode shapes. The alters are setup so that the 
final mass and stiffness matrices generated in the large displacement analysis 
are used. The advantages of utilizing the combined analysis capability are 
reduced CPU, faster turnaround time, and reduced quantities of output listings. 
In general, the CPU time can be reduced by one half. 
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APPENDIX A 


MSC/NASTRAN SOLUTION 64 INPUT DATA 


$ CRAY JOB CONTROL LANGUAGE (JCL) 


JOB , JN=BLADE64 , T=1 00 , MFL=450000 . 

ACCOUNT, AC=XXXXX t APW= XXXXX. 

NASTRAN , I D=XXXXX , DBASE=BLADE , FL=40000Q , STORE=TPOOL . 

/EOF 

$ NASTRAN EXECUTIVE CONTROL DECK 

ID BLADE, BLADE 
APP DISP 
SOL 64 
TIME TOO 
$ 

$ DMAP ALTER FOR COMPUTING CENTRIFUGAL SOFTENING MATRIX 
$ 

ALTER 77 

PARAMR //MPY/V , N , OMEGA/V , Y , R PM/ . 1 047 1 98 $ 

PARAMR //MPY/V , N , OMEGASQ/OMEGA/OMEGA $ 

PARAMR //C0MPLEX//0MEGASQ/0 . 0/V , N , ALPHAC $ 

MATMOD K1 , , , , ,/Kl B0B,/5/LUSET $ 

ADD K1 B0B,/K1 OM/ALPHAC $ 

MPYAD K1 OM,MJ J , /KS0FT/0/-1 //I /6 $ 

DBSTORE KSOFT, //MODEL/O/ $ 

$ 

$ DMAP ALTER FOR ADDING SOFTENING MATRIX TO STIFFNESS MATRIX 
$ 

ALTER 188 $ 

DBFETCH /HKSOFT, , , , /MODEL/O/1 $ 

ADD KJJ ,HKSOFT/KTOT// $ 

EQUIV KTOT,KJJ/ALWAYS $ 

ENDALTER 

CEND 
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$ NASTRAN CASE CONTROL DECK 


TITLE = SWEPT, CANTILEVER BLADE 
SUBTITLE = LARGE DISPLACEMENT ANALYSIS AT 8350 RPM 
LOAD = 50 
SUBCASE 1 

LABEL = LINEAR ANALYSIS 
SUBCASE 2 

LABEL = DIFFERENTIAL STIFFNESS 
SUBCASE 3 

LABEL = SUBCASE 3 
SUBCASE 4 

LABEL = SUBCASE 4 
SUBCASE 5 

LABEL = SUBCASE 5 
SUBCASE 6 

LABEL = SUBCASE 6 
DISP = ALL 
OLOAO = ALL 
SPCFO = ALL 
SUBCASE 7 

LABEL = SUBCASE 8 
DISP = ALL 
OLOAO = ALL 
SPCFO = ALL 
SUBCASE 8 

LABEL = SUBCASE 10 
DISP = ALL 
OLOAD = ALL 
SPCFO - ALL 
BEGIN BULK 


$ NASTRAN BULK DATA DECK 


C0RD2R 33 0 0.0 0.0 0.0 0.0 0.0 

1 .0 +2 

+2 -.866 0.50 0.0 

CQUAD4 1 10 1 24 3 

CQUAD4 2 10 3 46 5 

CQUAD4 3 10 5 68 7 

DMI K1 0 6 1 1 6 6 

DMI K1 2 2 1.0 

DMI K1 3 3 1.0 
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123456 

123456 


GRID 

1 

33 

-1.6 

GRID 

2 

33 

0.4 

GRID 

3 

33 

-1 .4 

GRID 

4 

33 

0.6 

GRID 

5 

33 

-1.2 

GRID 

6 

33 

0.8 

GRID 

7 

33 

-1.0 

GRID 

8 

33 

1.0 

PSHELL 

10 

10 

0.10 

MAT1 

10 

1 .0+7 


PARAM 

MAXRATI02.0+15 


PARAM 

RPM 

8350. 


PARAM 

K6R0T 

1.0 


RFORCE 

50 

0 


ENDDATA 




/EOF 





0.0 

5.9 

0.0 

6.1 

0.0 

3.9 

0.0 

4.1 

0.0 

1 .9 

0.0 

2.1 

0.0 

-.1 

0.0 

.1 

10 

0.30 

2.6 


139.2 1 .0 0.0 0.0 



APPENDIX 8 


MSC/NASTRAN SOLUTION 63 INPUT DATA 


$ CRAY JOB CONTROL LANGUAGE (JCL) 


JOB, JN=BLADE64,T=100,MFL=450000. 

ACCOUNT, AC=XXXXX, APW= XXXXX. 

NASTRAN , I D=XXXXX , DBASE=BLADE , ACCESS=TPOOL . 
/EOF 

$ NASTRAN EXECUTIVE CONTROL DECK 

ID BLADE, BLADE 
APP DISP 
SOL 63 
TIME 100 
CEND 


$ NASTRAN CASE CONTROL DECK 


TITLE = SWEPT, CANTILEVER BLADE 
SUBTITLE = FREQUENCY ANALYSIS AT 8350 RPM 
SET 1 = 0 
SEKR = 1 
SEMR = 1 
METHOD = 10 
DISP = ALL 
BEGIN BULK 


$ NASTRAN BULK DATA DECK 


C0RD2R 

33 

0 

0.0 

0.0 

0.0 

0.0 0 

1 .0 

+2 






+2 

-.866 

0.50 

0.0 




CQUAD4 

1 

10 

1 

2 

4 

3 

CQUAD4 

2 

10 

3 

4 

6 

5 

CQUAD4 

3 

10 

5 

6 

8 

7 

DM I 

K1 

0 

6 

1 

1 

6 

DMI 

K1 

2 

2 

1 .0 



DM I 

K1 

3 

3 

1 .0 




n 
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GRID 
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APPENDIX C 


CENTRIFUGAL SOFTENING 

The equation of motion for a lumped mass, In a rotating reference frame 
can be written ujtlng (ref 9): 

M (u> + 2 {Q} X {u} + {Q} x {Q} x [{R} + { u > ] = {F} (V) 

In this equation, M Is the mass, {U}, { u } , and (u } are the acceleration, 
velocity, and displacement vectors of the mass, {Q} Is the rotation vector, 

{ R } + { u > Is the position vector from the mass to the rotational axis, and 
{ F > Is the vector of applied forces. 


If the rotation Is about the Z axis then {Q} =01 + Oj + Ok and 
equation (1) can be written In matrix form as: 
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( 2 ) 


The third term on the left hand side of this equation contains the centri- 
fugal softening terms ( -fii 2 m) In addition to the elastic stiffness terms. From 
equation (1) It can be seen that the softening terms arise from the displace- 
ment u of the mass about It's Initial position. Since {u} normally Is much 
smaller than {R}, It Is tempting to neglect the softening terms, however, Q 
can be large, thus, for stability, large displacement, and normal modes analyses 
the softening terms are significant and cannot be neglected. 

An alternate way of deriving the softening terms x Is to consider the 
spring mass system subject to a centrifugal force £ 2 RM as shown In figure 6. 
The equilibrium equation for this system Is written: 

ku = Q 2 RM ( 3 ) 

where R Is. the distance from the mass to the stationary coordinate. Sub- 
stituting R = R + u Into the previous equation, the following Is obtained: 

ku = Q 2 (R + u)M 

( 4 ) 


(k - ft 2 M)u = ft 2 RM 

This approach produces the same softening terms found previously in 
equation (2). From this formulation It Is clear that the softening terms come 
from the fact that the centrifugal force is a function of both the Initial 
position R and the displacement u. Furthermore, the part that Is a function 
of u can be expressed as an Increase In the centrifugal force or as a nega- 
tive stiffness or softening term. 

When performing the large displacement analysis, NASTRAN will automati- 
cally consider the Increment In load by updating the centrifugal loads based 
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on the deformed position of the blade ( PARAM,SKPLOAD) . Although this approach 
will work for the large displacement analysis, the softening terms must still 
be Inserted into the stiffness matrix when a subsequent normal modes analysis 
Is to be performed. (In the sample problem both load updating and softening 
terms are used.) The softening terms must be Included In the stiffness matrix 
that Is transferred from the Solution 64 to the Solution 63 analysis, or the 
frequencies will be computed Incorrectly. 

To demonstrate the effect that the softening terms have on the natural 
frequencies, a steel plate (6" by 2“ by 0.10") was analyzed (fig. 7). In this 
figure two sets of frequencies are plotted; one for the plate lying In the 
plane of rotation, and the other lying perpendicular to the plane of rotation. 
For the latter case the softening terms have a significant effect on the 
plate's first bending mode frequency. This Is understandable since both the 
softening terms and the bending mode motion are In the plane of rotation. For 
actual blades, which have more complexity to their geometries than the steel 
plate, the centrifugal softening will have some Influence on all of the modes 
and will therefore need to be Included In all of the normal modes analyses. 
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APPENDIX D 


COMBINED SOLUTION 64/63 DATA DECK 

The following cards are added to the NASTRAN data deck for the combined 
analysis . 

Executive Control Deck 
ALTER 77 $ 

PARAM //SUB/V , N , MSUBS/V , Y , NSUBS/1 $ 

ALTER 311 $ 

PARAM //EQ/V , N , JEIG/MSUBS/NSK1 P $ 

COND LBL1 3,3EIG $ 

ALTER 315 $ 

DBFETCH /DYNAMICS, MJ J , GM, , /MODEL/O/O $ 

DBFETCH /GKAA, , , ,/SOLID/O/l $ 

MATMOD MJJ , , , , , /NVEC , /l 2/S , N , NULLS/2 $ 

EQUIV MJ 3 ,MNX/MPCF1 $ 

COND LBLB3,MPCF1 $ 

MCE2 GUSET,GM,MJ3 , „ ,/MNX, „ , $ 

LABEL LBLB3 $ 

COND LBLB1 ,NULLS $ 

SCE1 GUSET,MNX, , ,/MXBOB, , , , , $ 

JUMP LBLB2 $ 

LABEL LBLB1 $ 

EQUIV MNX,MXBOB/ALWAYS $ 

LABEL LBLB2 $ 

DPD DYNAMICS, GPLS»SILS, USET, SLT, /GPLD, SILD, USETD, , , , , , ,EED, 

EQDYN/LUSET/V,N,LUSETD/V,N,NOTFL/V,N,NODLT/V,N,NOPSDL/ 

V, N, NOFR L/ V, N,NONLFT/V,N,NOTRL/V,N, NOE ED/C, N,0/V,N,N0UE $ 

READ 6KAA,MXB0B, , ,EEO,GUSET,CASECC/LAMA, VECTOR, 

MI , OEIGS/MODES/S, N , NEIGS/NSUBS $ 

OFP LAMA,OEIGS// $ 

COND FIN,NEIGS $ 

SDR1 GUSET,, VECTOR,,,, GM, , , ,/UGV, ,QG/1 /REIG $ 

SDR2 CASECC , FCSTMS, FMPT, FDIT, FEQEXINS, , FETT, , FBGPDT, LAMA,QG, 

UGV , EST, XYCDB/OPG1 ,OQGl ,OUGVl ,OESl ,QEF2, PUGV/REIGEN/-1 $ 

OFP OUGV1 ,OPGl ,OQGl ,OEF2,OESl// $ 

LABEL FIN $ 

PRTPARM ////l $ 

ENDALTER 

Case Control Deck 

An additional subcase Is added to the subcases required for the normal 
application of solution 64. This subcase contains the METHOD ID of the EIGR 
card In the BULK DATA deck. 

SUBCASE XX $ LAST SUBCASE 
LABEL = EIGENVALUE EXTRACTION 
METHOD = 22 
OISP = ALL 
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Bulk Data Deck 

In addition to the EIGR card a new parameter NSUBS must be defined as the 
total number of subcases In the CASE CONTROL deck. 

PARAM, NSUBS, X $$ TOTAL NUMBER OF SUBCASES 
EIGR 22 
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FIGURE I. - TYPICAL FLEXIBLE BLADE FINITE ELEMENT MODEL. 
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FIGURE 2. 
BLADE. 
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BLADE ROTATIONAL SPEED, RPM 
FIGURE 3. - TYPICAL FLEXIBLE BLADE CAMPBELL DIAGRAM. 
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FIGURE A. - DEMONSTRATION PROBLEM FINITE ELEMENT MODEL. 



ROTATED BLADE DEFINED IN 
x',y\Z' COORDINATE SYSTEM 7 




FIGURE 5. - ALTERING ANGLE OF ATTACK USING ROTATED COORDINATE FIGURE 6 . - SPRING MASS SYSTEM UNDER CENTRIFUGAL LOADING. 

SYSTEM. 



FIGURE 7. - EFFECT OF CENTRIFUGAL SOFTENING ON BLADE 
FREQUENCIES (FIRST BENDING FREQUENCIES). 
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